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LONG-TERM  GOAL 

To  obtain  range-dependent  model  for  the  sound  speed  in  the  water  eolumn  and  sediment  properties 
from  the  dispersion  relation  of  the  propagating  modes. 

OBJECTIVES 

The  present  study  is  limited  to  the  estimation  of  range-dependent  sediment  properties  from  dispersion 
data. 


APPROACH 


The  acoustic  field  in  range  independent  shallow  water  can  be  represented  as  the  sum  of  contributions 
from  a  set  of  normal  modes  and  is  given  by  the  expression  given  below: 

P  {r  ,  z  ^  ,  z  ,co  )  =  in  Y,  4>  „  (  2  .  )4»  „  (  ^  ^  ^  ^  )  (1) 

n  =  1 

In  the  above  expression,  P  is  the  acoustic  pressure,  O)  is  the  frequeney,  Hq^  is  the  Hankel  function  of 
zero  order  and  first  kind,  N  is  the  number  of  modes,  (|)„  is  the  mode  function  of  the  nth  mode,  is  the 
eigenvalue  of  the  nth  mode,  and  z,.  are  the  source  and  reeeiver  depths,  and  r  is  the  range. 

Using  the  asymptotic  approximation  for  the  Hankel  function  the  above  equation  takes  the  form 


exp(  ikj  ) 
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In  a  range-dependent  ease  where  the  adiabatie  approximation  is  valid,  the  pressure  field  is  given  by  the 
expression 
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The  modal  phase  is  therefore  given  by  the  expression 


j  k  „  (  s  )  ds 
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Taking  the  differenee  on  both  sides  we  have 

r 
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Let  us  eonsider  the  situation  where  the  only  ehange  in  the  environment  is  in  its  eompressional  wave 
speed.  Then  at  any  given  range  loeation  s,  the  ehange  in  eompressional  wave  speed  is  related  to  the 
ehange  in  the  eigenvalue  by  the  relation  (Raj an,  1987) 
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where  Ae(s,z)  is  the  change  in  the  eompressional  wave  speed  at  range  s  and  depth  z,  c(s,z)  and  p(s,z) 
are  the  background  eompressional  wave  speed  and  density  respectively,  and  (])„  (s,z)  is  the  mode 
function  at  the  location. 


Inserting  this  expression  for  Ak^  in  equation  (5),  we  obtain  the  equation 
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The  travel  time  for  the  nth  mode  is  given  by 
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The  double  integral  can  be  changed  into  a  double  sum  given  below: 
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where  p  and  q  represent  the  range  and  depth  grids. 


Similar  expressions  relating  the  arrival  time  to  the  ehanges  in  the  eompressional  wave  speed  ean  be 
obtained  for  different  frequeneies.  This  set  of  equations  ean  be  redueed  to  a  matrix  equation  and  the 
equation  solved  to  determine  the  quantity  Ae  (s,z). 

The  experimental  set  up  for  this  approaeh  is  to  use  a  broadband  souree  for  transmitting  the  signal.  A 
reeeiver  plaeed  at  some  distanee  from  the  souree  aequires  the  signal.  By  performing  a  short  time 
Fourier  transform  or  wavelet  analysis  of  the  reeeived  signal,  the  dispersion  relationship  is  obtained. 
This  dispersion  relation  yields  the  arrival  time  of  eaeh  mode  at  different  frequeneies.  This  forms  the 
data  veetor.  The  elements  of  the  quantity  A  ( Sp ,  Zq)  are  eomputed  for  eaeh  mode  at  the  different 
frequeneies  for  the  initial  model  for  the  environment.  Having  eomputed  these  quantities,  the  matrix 
equation  is  solved  to  find  the  eorreetions  to  the  environmental  parameter.  As  the  problem  is  ill 
eonditioned,  regularization  methods  will  need  to  be  adopted  to  solve  the  inverse  problem. 

It  is  proposed  to  test  this  approaeh  using  data  to  be  aequired  in  two  field  experiments  seheduled  to  take 
plaee  in  the  years  2000  and  2001 .  The  first  is  an  experiment  sponsored  jointly  by  ONR  and  Japan 
Defenee  Ageney.  The  experiment  will  be  eondueted  in  the  general  area  of  the  Hudson  Canyon.  As 
part  of  this  experimental  program.  Woods  Hole  Oeeanographie  Institution  will  deploy  a  number  of 
drifting  sonobuoys  in  an  experimental  seenario  similar  to  experiment  that  have  been  performed  by 
WHOI  to  study  the  evolution  of  modes  in  the  shallow  water  environment.  Normally  this  experiment  is 
eondueted  using  a  souree,  whieh  puts  out  a  number  of  pure  tones.  A  broadband  eomponent  has  been 
added  to  this  experiment.  Two  types  of  broadband  signal  will  be  transmitted;  a  signal  that  is  a  sum  of 
CW  pulses  at  different  frequeneies  and  a  FM  sweep.  These  will  be  transmitted  when  the  sonobuoys 
are  at  range  of  10  -  15  km  from  the  souree.  The  signal  reeeived  by  the  sonobuoys  will  then  be 
analyzed  to  determine  the  dispersion  relation  for  the  modes.  We  propose  to  work  elosely  with  Dr.  G. 
V.  Frisk  (PI  for  the  WHOI  experimental  effort)  to  design  the  broadband  eomponent  of  the  experiment 
and  analysis  of  the  data.  The  seeond  field  experiment  is  the  East  China  Sea  experiment  seheduled  to 
take  plaee  in  the  year  2001.  During  this  experiment.  Dr.  James  Miller  from  the  University  of  Rhode 
Island  will  obtain  modal  dispersion  data  from  aeoustie  data  aequired  on  a  vertieal  array.  Explosives 
will  be  used  as  the  broadband  source  for  this  experiment.  We  propose  to  work  with  Dr.  Miller  and 
analyze  the  data  to  determine  the  range-dependent  bottom  properties. 

WORK  COMPLETED 

Two  types  of  broadband  signals  were  designed  for  the  US-Japan  experiment.  The  first  one  is  a  sum  of 
CW  pulses  and  other  is  an  EM  sweep.  The  dispersion  curves  for  the  two  types  of  signals  are  shown  in 
Eigures  1  and  2.  Both  these  types  of  signals  will  be  generated  using  J15  source.  In  the  experimental 
configuration,  the  source  will  be  stationary.  The  receiver  carried  by  the  sonobuoy  and  its  drift  during 
the  time  it  receives  the  broadband  signal  will  be  very  small  and  can  be  neglected. 

In  the  proposed  experiments  the  propagation  will  take  place  over  fairly  large  distances.  This  will 
require  the  range  dependent  inversion  to  include  inversion  for  the  water  column  sound  speed  profile 
and  the  bottom  depth.  The  linear  inversion  approach  outlined  above  can  be  modified  to  take  into 
account  changes  in  the  water  column  sound  speed  profile  and  bottom  depth  with  range. 


Source  signal  3  (ones  (30  Hz  50  Hz.  100  H2. 150  Ha  200  Hz) 


Figure  1.  The  transmitted  signal  consisting  of  5  CW pulses  and  the  dispersion  curve. 
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Figure  2:  Dispersion  curve  for  FM  sweep. 


However,  simultaneous  inversion  of  all  these  parameters  poses  a  ehallenge.  First  we  need  to  verify 
whether  the  problem  is  linear.  Simulations  eondueted  indicate  that  changes  in  modal  group  speed  with 
changes  in  bottom  depth  is  highly  nonlinear  where  as  the  changes  in  the  group  speed  with  changes  in 
bottom  compressional  wave  speed  is  only  weakly  non-linear.  The  non-linearity  of  the  change  in  group 
speed  with  changes  in  depth  is  weakest  for  the  lowest  mode.  This  analysis  show  that  it  is  necessary  to 
conduct  preliminary  bottom  surveys  to  get  good  initial  estimates  of  bottom  depth.  Further  in  order  to 
obtain  estimates  of  all  the  range  dependent  parameters  a  hybrid  scheme  that  uses  non-linear  approach 
and  a  linear  method  may  be  required.  Currently  inversion  with  simulated  data  is  in  progress. 

IMPACT/APPLICATIONS 

The  approach  outlined  above  permits  the  estimation  of  range-dependent  properties  of  the  environment. 
This  coupled  with  the  use  of  drifting  sonobuoys  can  be  used  determine  environment  parameters  in  real 
time  and  will  prove  an  excellent  asset  for  adaptive  sonar  design  and  for  rapid  environment  assessment. 
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